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Abstract 

The lepton-number violating process e^p — > i'(,{T^)l^l'^X mediated by Majorana 
neutrinos is studied for the HERA collider for (//') = (er), (/^r), (///i) and (rr). 
Only the muonic decay of the r is considered. The direct limit on the effective muon 
Majorana mass, {m^^) is improved significantly to 4.0 • 10^ GeV and for the first 
time direct limits on the analogous effective masses connected with the tau sector 
are given, namely 4.2-103 GeV for (mer), 4.4- 10^ GeV for {m^r) and 2.0-10^ GeV for 
i^mT-T)- We find that a more general analysis for an upgraded HERA could improve 
this values by a factor of up to 40, yet still being orders of magnitude worse than 
indirect limits. 
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1 Introduction 



Since there is growing experimental evidence of nonzero neutrino masses using neutrino 
oscillations as explanations 0] , an additional fundamental question still to be solved is the 
character of the neutrinos, i. e. are they Dirac or Majorana particles? From the theoretical 
side the latter case is favored since Majorana particles pop out of almost every GUT |^ 
and are also the product of the attractive see-saw-mechanism The most important 
tool to answer this question is the detection of lepton-number violation in the neutrino 
sector. The most effort of theoretical and experimental work has been put in neutrinoless 
double beta decay (Oz//3/3), resulting in an upper limit on the effective Majorana mass 
(rriee) = I U^j^rrimV^l of about 0.2 eV 0, where the mass eigenvalues, = ihl 

the relative CP parities and f/em the mixing matrix elements. In general, there is a 3 x 3 
matrix of effective Majorana masses, the elements being 

(m„^) = \ j2UamU(3mmmV^ with a, (3 = e, /x, r. (1) 

For sake of simplicity we assume that the elements Uam are real and skip also t]^ ■ Only few 
direct information on elements other than (mge) is available: muon-positron conversion in 
sulfur gives a limit on (m^e)~ 0.4 (1.9) GeV, when the final state proton pairs are in a spin 
singlet (triplet) state, respectively. This limit is obtained when comparing the theoretical 
value from with the PDG branching ratio limit and using the fact that the matrix 
element of the process is proportional to |(m^e)P- In a recent paper |^ we considered 
the reaction u^N — >• fi~ fj,'^ fi'^ X , mediated by Majorana neutrinos, and deduced a limit 
of (m^^) ~ 10"^ GeV, improving the previous bound by one order of magnitude. To 
our knowledge, there are no direct limits on other elements of (m^^). Note that we are 
considering direct limits, i. e. using processes sensitive on the respective quantity. Indirect 
bounds, obtained from oscillation experiments and unitarity of the mixing matrix will of 
course be far more restrictive. 
In this paper we will study the process 

e^p IH'^X, with (//') = (er), (/ir), (/i/i) and (rr) (2) 

for the case of the HERA collider. We will focus mainly on the e^p mode, but the qualita- 
tive conclusions we draw remain of course valid for e~p-collisions as well. We demand the 
tau^ to decay in muons to take advantage of the like-sign lepton signature, which is more 
unique for muons than electrons. A further effect is to have a small number of like-sign 
(e/x) or (ee) events, which might be more background dominated. The relevant diagram 
for process (Q) is shown in Fig. |1|. 

It is evident that such a process has a spectacular signature with large missing transverse 
momentum {pj) and two like-sign leptons, isolated from the hadronic remnants. Direct 



production of heavy Majoranas at HERA has been studied before |jTO[ with the process 
e~p — >■ XN — i> XW^P Xul'^P resulting in two leptons with different charge. In 



^We shall use the term electron, muon or tau for both, particle and antiparticle. 
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Figure 1: Feynman diagram for the process e^p l^l'^X. q denotes the four- 

momentum of the propagating Majorana neutrino. Note that there is a crossed term 
and for / 7^ /' two possible assignments for the corresponding lepton vertex exist. 



contrast to direct production the process discussed in the present paper deserves some 
attention because of its unique signature. 



2 Analysis 

We use HERA kinematics for electron and proton energies = 27.5 GeV, Ep = 820 GeV 
and the GRV 98 |TT| set of parton distributions. To mimic the experimental situation 



the following kinematical cuts are applied: \rii^x\ < 2.0 and > 10 GeV. Here 



{J2Pxy + (Sp?^)^ is the total missing transverse momentum, with the sum going over all 
neutrinos, i. e. one neutrino for the (/i/i), three for the (er) and (fir) and five for the (rr) 
case, rji^x = — lntan(6';^x/2) is the pseudorapidity of the charged lepton / and the hadronic 
final state X respectively, with 6 the polar angle in a system where the 2;-axis is parallel 
to the proton direction. In addition, we want the charged leptons to be isolated from the 



hadrons and demand AR = \J{4>i — 4>xY + {Vi ~ VxY > 0-5 with being the azimuthal 
angle. 
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The exact calculation of the diagram and some general features of the resulting cross section 
are described in P]. For the problem at hand we additionally folded in the three-body 
decay of the tau leptons. When considering heavy neutrinos, one has to note the mixing of 
the usual standard model (SM) leptons with these hypothetical particles. The Lagrangian 
for the lepton I coupling to neutrino mass eigenstates Um is: 



V m 

This leads to a dependence of the cross sections of the form: 



(3) 



da{ep — > Uell'X) oc 



nir. 



I'm ' 



mt 



(4) 



with q being the four-momentum of the propagating Majorana neutrino. 

ViVrn and found a limit of |f//mP < 
80 GeV and I = e, fi and r. For larger masses analyses of 



The DELPHI collaboration [|I2| examined the mode Z 
2 ■ 10~^ for masses up to 



neutrino-nucleon scattering and other processes yielded 

^ \Uem\^ < 6.6 ■ 10-^ \U^.m\^ < 6.0 ■ 10"=^ and 1^™!' < 1-8 ■ lO'l 



(5) 



3 Results and Discussion 

In Fig. ^ the total cross section as a function of one mass eigenvalue rrim is shown for all 
combinations of final state charged leptons, without considering the above mentioned Uim 
limits (i. e. setting = !)• Our condition that the tau decays into a muon is included. 

As can be seen, there is a maximum at about 70 GeV, which has purely kinematical 
reasons, see [§]. This means that we can assume one eigenvalue dominating the sum 
S ^'ini'f ~ ^m)~^- For small masses the cross section rises with and for higher masses 
it falls with as can be understood from the two extreme limits of m^(g^ — m^)~^. 
The masses of the final state leptons have almost no effect (less than 5 %), so that the 
only numerical difference comes from the branching ratios (BR), the Uim limits from Eq. 
(^) and the factor 2 for the (er) and (/ir) cases. The latter comes from the two possible 
assignments the two leptons have, when emitted from the Majorana vertex. 

These limits combined with the branching ratio (BR(r — > iii^^fT) = 0.1732 0) lead to 
a maximal cross section of about 10~^^ b in the (/z/z) case for a neutrino mass of 80 GeV. 
It is about 12 orders of magnitude smaller than the SM charged current (CC) e'^p cross 
section of acc(e^p, > 200 GeV^) ^ 30 pb (we checked that the Q"^ condition is not 
significantly violated for the cuts we applied). Nevertheless the above cross section is some 
orders of magnitude closer to the relevant SM CC process than most other exotic Majorana 



neutrino induced AL = 2 processes such as — *• vr JTSf or n /i """-conversion via 

muon-capture in ^"^Ti [|TB|. These have ratios with respect to the relevant SM CC process 



of at most 10 For u^N — > fi fi~^fi~^X a ratio of 10 ^"^ [§] for a 500 GeV neutrino beam 
is achieved, so the process described here results in an improvement of another 5 orders of 
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Figure 2: Total cross section for e^p V^l'^l'^X as a function of one eigenvalue nim- No 
limits on Uim are applied, the branching ratio for taus into muons is included. The (er) 
and (/ir) cases are indistinguishable in this plot. 

magnitude. 

The e~p mode gives (in contrast to the normal CC process) a cross section smaller than 
for e^p by a factor of about 2 (1.9 for small masses, 2.5 for masses higher than 10^ GeV), 
so that the ratio for this mode is a factor 4 worse. The other cases, like (rr), have ratios 
with respect to the SM CC process maximally one order of magnitude smaller than the 
(/i/i) case. 

As an example for differential cross sections we plot in Fig. |^ for = 80 GeV the dis- 
tribution of the missing transverse momentum for the (/i/i), (/ir) and (rr) case. Note 
that all these cases have two like-sign muons in the final state. The mean values are 
(j^) ~ 28.3, 37.0 and 37.3 GeV, respectively. The shape is different for each case, de- 
spite the similar mean value of missing transverse momentum. To distinguish, say, (/ir) 
from (rr) events, one should consider other distributions, e. g. the invariant mass of the 
two muons, m(2/i), as displayed in Fig. ^. Whether a muon comes directly from the 
Majorana vertex or from the tau-decay makes its energy and momenta fraction of the 
total available energy smaller and changes its invariant mass. Here the mean values are 
(m(2/i)) ~ 66.8, 28.1 and 16.1 GeV, for (/i/i), (/ir) and (rr), respectively. Unfortunately, 
this procedure requires high statistics. On an event-by-event analysis detailed kinematic 
reconstruction as well as angular distributions might be more useful. 
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Figure 3: Distribution of the total missing transverse momentum ^ for the (/xr), (/x/x) and 
(rr) case in the reaction e+p — > V^l'^V'^X for rrim — 80 GeV. In order to have the curves 
in the same order of magnitude we did not include the tau branching ratio. 
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Figure 4: Invariant mass distribution for the two like-sign muons given for the same 
parameters as in the previous figure. 
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Because of its lepton-number violating character the discussed process should be back- 
ground free from standard model processes. On the other hand, misidentification of the 
muon charge might create some like-sign dimuon events. Processes producing opposite 
sign muon pairs are pair production of heavy quarks, photon-photon interactions, Z-boson 
production, Drell-Yan pairs from resolved photon-proton interactions, pion punchthrough 
associated with a single muon event and beam related background. Practically all these 
processes can be eliminated by kinematical arguments. Therefore the identification of two 
like-sign, isolated muons with large pt in addition to a large missing transverse momentum 
should indeed be an outstanding signature. For the (er) channel, which has a like-sign 
efj, signature, single W production is a severe background. However, sophisticated 
kinematical arguments as those given in [|l^ for the case of direct Majorana production 
might be also applicable in the case discussed here. Because of this we shall assume zero 
background, leaving a more detailed analysis for further studies. In case of observation 
such a detailed analysis has to be done anyway in order to rule out any standard model 
process as the ones described above. 

In order to get bounds for {niap) we assume that the cross sections displayed in Fig. ^ 
are proportional to |(?tIq,^)P and take the luminosities used in searches for isolated lepton 
events with missing transverse momentum, i. e. £e+ = 36.5 pb~^ (HI, [^) and £e+ = 47.7 
pb~^ (ZEUS, [|19|). We take the average of 42.1 pb~^ and get values in the range of 10^ 



to 10^ GeV, thereby improving the (rn^^) limit with respect to [§] significantly and giving 
for the first time direct limits on (mar)- Combining all limits, ignoring possible phases in 
the elements Uam (therefore getting a symmetrical matrix {rriais)) as well as skipping the 
intrinsic CP parities, the following direct bounds for the effective Majorana mass matrix 
exist: 



' {niee) {rrie^,) {rrier) \ 



/ 2 • 10-1° 0.4(1.9) 4.2 ■ 10^ \ 



< 



V 



4.0-103 4.4-103 

14 



GeV. (6) 



2.0 - 10^ / 



A spread over 14 orders of magnitude can be seen. We state again that these are direct 
limits and the elements other than (rriee) should not be confused with their real values. 
As is evident and not surprising, the bound coming from {)v(3(3 is by far the best limit for 
an effective Majorana neutrino mass. One might argue that FCNC processes like r — /i7 
place severe bounds on this effective masses. Applying the BR from |^ to the measured 
limits from 0, ^ gives bounds for = 



of the order 1 to a few 10 GeV. 



Without specifying to a special mass and mixing scheme it is rather difficult to compare 
{rriais) with m^. Since there is no commonly accepted scheme around, we believe that 
numbers derived from experiments are necessary. Another point is that in principle one 
could derive the remaining elements of {map) from the limits on (mg^) and (mee). Here 
the same argument holds. One should say that if the value for (mee) is fixed, all other 
elements of {niap) should be in the same order of magnitude, therefore at most a few eV. 
Our matrix (^ is thus far from being physically realized. 

The factor of about 3 the (rr) limit is worse is due to our condition that only the tau decay 
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into muons is considered, which could be skipped in a more general analysis including more 
decay channels. Note that with our assumption a oc {rriaiiY the bound is proportional to 

^JxfoC. Therefore a general treatment of all possible tau decay channels would bring a 
factor of about 5.8 for (rr) and 2.4 for the channels involving only one tau. Furthermore, 
an optimistic luminosity value of 1 fb~^ would bring another factor of about 5. 
An upgraded HERA with — 1020 GeV, E^. = 33.5 GeV rises our cross sections for the 
small mass regime by about 60 %, so that with an integrated luminosity of 1 fb~^ and in 
consideration of all tau channels our bounds could thus be lowered by a factor of 40. 

4 Conclusions 

We have studied the Majorana neutrino induced process e^p — > W l^l'^X at HERA and 
deduced for the first time direct bounds on all effective Majorana masses other than the 
one measured in 01/(3(3. A way to distinguish signal events from each other as well as from 
background is discussed. We propose a search for two like-sign muons in the final state 
combined with large missing transverse momentum. The cross sections are typically 12 to 
13 orders of magnitude smaller than SM CC processes, which has to be compared with 
related rare meson decays or conversion on nuclei, which give ratios of at most 10~^°. 

We improved the direct limit on (m^^) significantly and gave for the first time direct limits 
on (mer), (TTifxr) ^ud (771,-^). Howcvcr, these are direct bounds, which will be orders of 
magnitude worse than ones derived indirectly. 
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